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Abstract 

Results of a calculation of the c-axis infrared conductivity o~ c for a d x 2_ y 2- 
wave superconductor which include both elastic impurity and inelastic spin- 
fluctuation scattering are presented and compared with the afr-plane conduc- 
tivity Gab in the same model. In this model, the interlayer c-axis coupling is 
taken to be weak and diffusive. While in clean systems, inelastic scattering 
leads to a peak at u> = 4Ao in a a b for T < T c , it has little effect on the 
corresponding a c , which exhibits structure only at u) — 2Aq and is directly 
related to the single-particle density of states N{u). The c-axis penetration 
depth A c in the same model is predicted to vary as T 3 at low temperatures 
in clean samples. We discuss recent optical experiments on the cuprates and 
compare with these predictions. 
PACS: 74.25.Nf, 74.72.-h, 72.10.Di, 72.10.Fk 
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I. INTRODUCTION 



In the superconducting state of optimally doped YBa 2 Cu 3 7 ^s, a gap-like depression 
develops in the frequency dependence of the c-axis conductivity for T less than T c . As em- 
phasized by Homes et al, EH this behavior is quite different from the in-plane ab conductivity 
observed in samples from the same source. We have recently developed a simple model for 
the ab conductivity which provided a qualitative fit to experimental data, and here we ex- 
tend this model to see how well it can describe the c-axis conductivity. The model consists of 
BCS quasi-particles with a d x 2_ y 2 gap and a lifetime determined by both elastic and inelas- 
tic scattering processes. The latter inelastic lifetime plays a key role in our model and was 
calculated from a spin-fluctuation interaction in which the quasi-particle d x 2_ y 2 gap entered 
explicitly in determining the spectral weight of the spin-fluctuations. The elastic scattering, 
which is important at low frequencies, was treated in the resonant impurity scattering limit. 
Here, in discussing the c-axis conductivity, we will treat the interlayer transport as weak 
and diffusive. Our conclusion from this analysis is that a model of a d x 2_ y 2 superconductor 
which includes inelastic as well as elastic scattering can provide a sensible fit of both the ab 
and the c-axis conductivities. In particular, it predicts the existence of absorption features 
at 4A in a a ^ and at 2A in a c , where A is the (i-wave gap maximum; such features are 
apparently observed in a variety of materials. To the extent that the interlayer coupling can 
be treated as weak and diffusive, the c-axis conductivity provides direct information on the 
single particle density of states in a similar manner to SIS tunneling. 



II. C-AXIS CONDUCTIVITY 

The nature of c-axis transport in the cuprates is still controversial, but it is widely 
believed to be weak and incoherent .H§ A simple model which describes this has H = H +H c , 
where H describes free particles with residual d x 2_ y 2 BCS pair interaction, and 

H c = Y,V l {c\ la c l2cT + h.c.) 1 (1) 
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where the V* may be thought of as randomly placed static impurities mediating hopping 
between planes 1 and 2. We emphasize, however, that the end result of our treatment will 
yield a description in terms of a single c-axis relaxation time which must be similar to any 
weak diffusive transport model. The interplane paramagnetic current operator for a given 
configuration of the random potential is then 

f^-ied^Viicl^-h.c). (2) 

i 

The response to an applied electromagnetic field A is now determined by expanding the 
kinetic energy of the system up to second order in A, yielding for the expectation value of 
the total current 

(j c (q,co)) = [d 2 e 2 (H c ) + A cc (q, cu)]A z (q, u), (3) 

where (jI)a = A CC A Z to leading order. 

We first examine the conductivity cr c (q = 0,u), assuming that V is weak, thereby ne- 
glecting multiple scattering vertex corrections as well as self-energy corrections, since one 
power of V is already contained in each of the c-axis current operators which enter A cc . 
After impurity averaging, we find 

a c (uo) = , (4) 

oo 

with 

lmA cc (uo) = 

4 tt d 2 e 2 J duo' £ |Vp„ p /| 2 (f(u + a/) - /(a/)) (5) 
pp' 

x (A(p, to + co')A(p', to') + B(p, to + co')B(p', to')) . 

Here the normal and anomalous spectral functions are defined as A{p,uo) = 
-^ImTr{G(k, oo + iO)}, B(p,oo) = -^ImTr{riG(k, uo + iO)}, where G is the matrix Green's 
function 

^ ^ = u2_C2_ A 2 ■ ( 6 ) 
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The interlayer scattering matrix element V p - p i is defined by 

Vi = Y, W - Rm) = E V p e^< n >- n ™\ (7) 

where m runs over all lattice sites and i labels an impurity site. For simplicity we take 
a cylindrical Fermi surface with £ p = (p 2 /2m — fi) and a d x 2_ y 2 gap A p = A (T) cos 2 0. 
The renormalized frequency uj includes the effects of scattering associated with the in-plane 
impurities and the lifetime effects due to dynamic spin fluctuations. In the resonant limit, 
the in-plane impurity scattering gives a renormalized quasi-particle energy u determined by 
the self-consistent solution of 

u = u- S imp (u;) = uj + iT/ ( ) . (8) 

Here T = rii/TrNo, with rii the impurity concentration and N the Fermi level density of 
states. 

We assume further that the inelastic scattering rate r i ~ 1 (c<j,T) arising from spin fluctua- 
tions can be calculated separately, and take as a model the one-antiparamagnon self-energy 
calculated in a 2D Hubbard model by Quinlan et al% The momentum dependence of the 
full l/rj n (p, u) has been neglected here in favor of a simpler effective relaxation rate for 
quasiparticles at the gap nodes, l/r in {uj) = 1 / T in (k noi ie, oj) as discussed in Refs. and |8| At 
low temperatures Tj" 1 varies as to 3 for u < 3 A , reflecting the usual uj 2 dependence of an 
electron-electron scattering process along with a third factor of to arising from the d x 2_ y 2 
density of states. For u > 3 A , 1/T in crosses over to the pseudo-linear behavior found in 
spin-fluctuation scattering calculations in the normal state. Although for concreteness we 
have chosen a particular model for l/r^, we note that the behavior in the superconducting 
state is characteristic of any model based on electron-electron scattering with quasiparticles 
which are well-defined at sufficiently low temperatures, in the presence of gap line nodes. 

Combining the elastic and inelastic contributions, we set 



in Eqs. (8) and (9). The resulting lifetime for in-plane scattering is plotted in Fig. 1 for 
several values of the impurity concentration. Note the Hubbard parameters have been chosen 
so as to give a value of l/r in (u; — > 0, T c ) = 2T C . 

As discussed by Graf et a/.,i setting |Vp_ p /| 2 equal to a constant corresponds to diffuse 
transmission, while taking an extreme forward scattering form such that p = p' (that is, 
= 0') corresponds to specular transmission. In the diffuse transmission limit, the second 
term in Eq. (5), arising from pair transfers, vanishes for a <i-wave gap. In the actual system, 



V pp i\ 2 is anisotropic, and here we model this using the separable form 



V pp ,\ 2 = \V\ 2 + \V'\ 2 cos 2 cos 2 0' 



(10) 



The interaction strength V can be adjusted to fit the normal state c-axis conductivity 
a cn = 2d 2 e 2 N /r c with r~ x = 2nf7rN \V\ 2 . Here n\ ~ 0(1) is formally the density of 
scattering sites which allows the planes to communicate. The second term V can be used 
to fit the c-axis penetration depth A c (0) (see below) 
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For a n ~ 150 (ficm)^ 1 and A c (0) ~ 10 4 A, we find \V /V\ 2 = 0(1). 

With the form (10) for |VJp_ p /| 2 , the expression Eq. (4) for a c (u) simplifies to 



(11) 
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x N(lu + lu')N(lu') + 
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N(u) = Re 



Co 2 - At 
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-Re K(A /cD) 
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(14) 



-Re{ — [K(A /i»-E(Ao/a»]}. 
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Here K and E are the complete elliptic integrals of the first and second kinds, respectively.0 
Results for a c (uj)/a cn versus to at different temperatures are shown in Fig. 2 under the 
assumption of isotropic diffusive scattering ( | X^'/V"| 2 = 0). 

The primary feature seen in the figure is the pseudogap of width 2Ao which opens up 
below T c . The low-frequency behavior can be estimated in this limit to be a c (u)/a cn ~ 
2o; 2 /3Aq. Figure 3 illustrates the effects of adding anisotropy to the interplane scattering, 
as well as intraplane scattering effects. Because of the form of Eq. (12) as a convolution of 
normal and anomalous densities of states, the shoulder feature at 2Ao(T) is quite robust. 

While strongly scattering in-plane impurities (r > 0) can alter the very low-frequency 
behavior, the 2A shoulder is essentially unaffected. The large intraplane relaxation rate 
near T c might be expected to change the behavior in an important way, but this is seen 
not to be the case, since the effect of a large 1/r is simply to smear out the density of 
states factors entering Eq. (12). Scattering anisotropy (IV'/V) 2 > 0) also does not lead to 
important qualitative differences from the results of Fig. 2, since the anomalous terms in 
the convolution give rise to a low-energy contribution to the conductivity of order (cu/Aq) 4 . 
The results in these limits are therefore similar to those obtained earlier by Graf et al.i In 
Fig. 4 we compare these results to the experimental data of Homes et al; here we have 
chosen a parameter set which appears to give a reasonable fit to the ab penetration depth 
and conductivity, but the c-axis conductivity is not terribly sensitive to this choice. 

Due to the weak behavior of the anomalous density of states at low temperatures, we 
note that within this model a measurement of the c-axis conductivity provides a nearly 
direct measurement of the in-plane density of states N(u). It will be interesting to compare 
N(u) obtained in this way to results obtained, e.g., from tunneling measurements. 

III. C-AXIS PENETRATION DEPTH 

The real part of the current-current correlation function may also be calculated, yielding 
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A cc (q ^0,u = 0) = -2e 2 d 2 <T^ £ IVp'P (15) 

<+>n pp'a 

x ^Tr[roG (T (p,^ n )r G (7 (p / ,^„)]. 
The diamagnetic term in Eq. (3) takes the form 

(H c ) = 2n t Tj2EK- P >\ 2 (16) 

<+>n pp'a 

x ^Tr[r 3 G (J (p,w n )r3G (7 (p , ,^ n )], 

which means that only the product of the anomalous Green's functions survives in this 
approximation. 

The remaining frequency sums and integrals may be evaluated in the diffusive limit 
treated above to find the London penetration depth, A C (T), as 

° -a cn \V'/V\ 2 J du tanh(^) Im[$H 2 ] (17) 



a cn |V7V| 2 A (T)[0.480 - ^(-r^)\ (18) 



4vrA c (T) 5 

~ n iT/'/T/i 2 A„mrn 4«n - - 

vr ^A (T)< 

where $ = |^-[K(Ao/£) — E(Aq/u))] is the same angular average appearing in Eq. (14). 
Note the last approximate equality is valid for T <C T c in the clean limit only. The expression 
(17) is equivalent to the result obtained by Radtke et a/.i in the absence of direct interlayer 
tunneling. In Fig. 5 we have plotted the temperature dependence of A c (0) 2 /A c (T) 2 for 
various values of the in-plane impurity scattering rate. As in the ab case,@ we expect a 
crossover to a 5X C ~ T 2 behavior at low temperatures when the system becomes sufficiently 
dirty. Although the normalized superfluid density (A c (0)/A c (T)) 2 is quite insensitive to the 
presence of planar impurities, as seen in Fig. 5, we note the peculiar result that in this model 
the penetration depth itself scales roughly with the c-axis conductivity, and should therefore 
be quite sensitive to disorder which affects c-axis transport. 

IV. COMPARISON OF AB AND C-AXIS DATA 



In constructing phenomenological theories of this kind, it is important to see if a consis- 
tent description of a wide class of experiments can be obtained, along with an understanding 
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of the effects of disorder. In this regard we would like to predict correlations among vari- 
ous quantities which should hold at least for all optimally doped materials, and possibly in 
the underdoped systems as well. The most important such correlation to emerge from the 
current analysis is the prediction of a clear shoulder in the c-axis conductivity at u — 2Ao, 
a broad maximum in the ab conductivity at 4A , and a shoulder-like feature in 1/T a f,(uj) 
at 3Ao- Since the gap maximum Ao is not known for all materials,^ one can search for 
correlations of such features in optical experiments. While such a method will not provide 
the most accurate measure of the gap size, it will serve to confirm the general picture of 
a d-wave order parameter with nodal quasiparticles scattered by resonant defects of some 
kind at low energies, and strong spin-fluctuation scattering at higher energies. In Fig. 6, we 
show er c , a a b, and l/r ab together to illustrate the correlations we expect. In Fig. 7, we show 
data on YBa^Cu^O^^^ YBa 2 Cu 3 OQ^, and Y Ba^Cu^O^, from recent work by Basov et a/.lll 
Here we would like to note some similarities between the results for the optimally doped 
material below its T c , and the behavior of the underdoped materials above their transition 
temperatures in the so-called pseudogap regime. The most distinct feature is the shoulder 
in a c (upper panel of Fig. 7), which, as we have discussed, occurs at 2A ~ 6 — 8T C for 
Y Ba<iCuzO%$§. In the underdoped compounds this shoulder, at what we will call "2A ", 
becomes visible above T c in the pseudogap regime. In the underdoped compounds, "2A " 
from this criterion is slightly reduced from its value in the optimally doped materials. The 
second panel in Fig. 7 shows o^iuS). In Ref. ||, the absence of a 2A structure and the broad 
rise at 4A observed in the afe-plane conductivity of YBa 2 Cu 3 O e ,g 5 below T c at 1000 cm" 1 
was attributed to a 4Ao inelastic absorption threshold smeared by (i-wave nodal quasiparti- 
cles. Now, as is well known, there is no gaplike structure in the ab conductivity (second panel 
of Fig. 7) in the doping range for which the pseudogap is observed in the c-axis conductivity; 
there is, however, a narrowing of the Drude peak and a broad rise at higher frequencies near 
"4Ao". It could be that the appearance of the "2Ao" in a c but not in a a b is similar to 
the behavior we have found for these quantities below T c in the model we have applied to 
YBa2Cu^OQ,g5, in which the crucial factor is the direct relevance of the density of states 
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N(w) to c-axis transport. Support for the correlation of the c-axis pseudogap and the den- 
sity of state is found in studies by Homes et a/.0 which found scaling of the pseudogap with 
the Knight shift in underdoped cuprates. Finally, as shown in the bottom panel of Fig. 7, 
that the in-plane relaxation rate for the underdoped materials shows a strong supression of 
the linear frequency dependence in the pseudogap regime even above T c . In all cases shown, 
the crossover to stronger frequency dependence at sufficiently low temperatures occurs at 
an energy consistent with the "3Ao" smeared threshold expected from the current model. 

We turn finally to the comparison of the c-axis penetration depth with its ah counterpart. 
The results are shown in Fig. 8 for the clean limit, and compared to the data of Hardy et 
a/.@ on YBa^Cu^O^^. It is clear that the result agrees qualitatively with the much flatter 
temperature dependence of the c-axis penetration depth found at low temperatures. The 
low temperature behavior is not yet well characterized, however. The discrepancies result, 
in the intermediate temperature range, from the crude isotropic band structure assumed, 
and in the range close to T c from critical fluctuations not accounted for in the current model. 

V. CONCLUSIONS 

We have argued that a simple theory of <i-wave superconductivity, together with impurity 
and spin fluctuation scattering, accounts for most of the qualitative aspects of both optical 
conductivity and London penetration depth measurements for currents along both the a and 
c axes, c-axis properties have been calculated under the further assumption that transport 
is primarily diffusive. The optical conductivities a c and a a b may then be correlated in this 
model, in that we expect a broad peak in a a b at u ~ 4A and a shoulder-like feature in 
cr c at uj ~ 2A . This behavior indeed appears to be reproduced qualitatively in several 
materials, including underdoped materials, where a behavior qualitatively similar to that 
observed in the optimally doped materials below T c is seen significantly above T c . If the 
phase coherence implicit in the BCS approach is ignored, the model considered here offers 
a natural interpretation of these results and an explanation of why the normal state c-axis 
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pseudogap does not manifest itself in planar transport properties. Further study is needed 
to put these speculations on a firmer basis. 

Finally, we have calculated and compared the c-axis London penetration depth in the 
same model, and shown that in clean systems it will vary as 5\ c ~ T 3 at low temperatures. 
This temperature dependence is much stronger than the ab result, 5\ a b ~ T in the same 
model, in qualitative agreement with experiments on optimally doped Y BaiCuzO-j-?,. 
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FIG. 1. Normalized intraplane inelastic relaxation rate Tj~ (o>)/Ao(0) vs. cj/Ao- Solid line: 
T = 0; dashed-dotted line: T = 0.8T C ; dashed line: T = T c . In all cases A (0)/T c = 3. 
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FIG. 2. Normalized c-axis conductivity a c (u)/a cn vs. w/Ao(0). Solid line: T = 0; 
dashed-dotted line: T = 0.5T C ; dotted line: T = 0.8T C ; dashed line: T = T c . Scattering anisotropy 
parameter |F'/y| 2 =0, planar impurity scattering rate T = 0, Aq(0)/T c = 3. 
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FIG. 3. Normalized c-axis conductivity a c (uj)/a cn vs. cj/Ao(0) for T = 0. a) |^-| 2 = with 
scattering (dashed line, V = 0.1T C and inelastic scattering included) and without scattering (solid 
line, 1/r = 0); b) 1/r = with |^| 2 = 1 (dashed line) and |^| 2 = (solid line). 
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FIG. 4. Normalized c-axis conductivity o c (uj)/a cn vs. u>/A (0), \V /V\ 2 = 0, T = 0.02T C , 
inelastic scattering included, T = T C (dashed line), T = 0.8T C (dashed-dotted line), and T = 0.1T C 
(solid line), for Aq/T c = 3. Data from Homes et al% on YBa2Cu^OQ,^, solid line: T=10K, dashed 
line: T=100K. u cn normalized to 100K data at 100 cm 
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FIG. 5. Inverse c-axis London penetration depth A c (0) 2 /A c (T) 2 vs. normalized temperature 
T/T c . Solid line: intraplane impurity scattering rate T/T c = 0; dashed-dotted line: T/T c = 0.05; 
dashed line: T/T c = 0.1. Inelastic scattering included. 
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FIG. 6. a) Normalized c-axis conductivity <j c {uj) / a cn vs. w/Ao(0); b) a6-plane conductiv- 
ity Pab{w)/<JQ vs. w/Ao(0) (Note ao = cr a b{u = 0) at T = T c ); c) quasiparticle relaxation rate 
r i ; 1 (w)/A vs. w/A (0). All panels show results for T/T c = 0.02, T = 0.1T C (solid line), 0.8T C 
(dashed-dotted line), and T c (dashed line). 
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from left to right. Dashed-dotted lines: 300K, dashed lines: T ~ T c , solid lines: 10K. The c-axis 
conductivity of single crystals is taken from Homes et al. and multiplied by a factor of 4. Lines 
indicate position of shoulder (2Ao) in a c , and twice this energy. (Lines not included in LANL 
archived preprint.) 
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FIG. 8. Inverse c-axis and afc-plane London penetration depth \ c ,ab(fy 2 1 ^c,ab(T) 2 vs. normalized 
temperature T/T c for intraplane impurity scattering rate T = 0, no inelastic scattering included. 
Data from Hardy et aZ0 on Y Ba2CuzO§$5. 
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